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We report on the evolution of tropospheric nitrogen dioxide (NO2) over Spain, focusing on the densely 
populated cities of Barcelona, Bilbao, Madrid, Sevilla and Valencia, during 17 years, from 1996 to 2012. This 
data series combines observations from in-situ air quality monitoring networks and the satellite-based 
instruments GOME and SCIAMACHY. The results in these five cities show a smooth decrease in the NO2 
concentrations of —2% per year in the period 1996-2008, due to the implementation of emissions control 
environmental legislation, and a more abrupt descend of —7% per year from 2008 to 2012 as a consequence 
of the economic recession. In the whole Spanish territory the NO2 levels have decreased by —22% from 1996 
to 2012. Statistical analysis of several economic indicators is used to investigate the different factors driving 
the NO2 concentration trends over Spain during the last two decades. 



Ozone, fine particulate matter (PM), nitrogen oxides (NOx = NO + NO2), and sulfur dioxide (SO2) are 
the most important air pollutants in the troposphere, which represent a serious threat for ecosystems and 
human health'. Although the sources of these pollutants in the European Union (EU) have considerably 
been reduced in the last ten years, the air quality is still far from the levels required by the EU legislation. The EU 
environmental policies requirements are currently being reviewed to be adapted to safer levels for human health''. 

NOx are important in tropospheric chemistry for several reasons, including i) their role in the photochemical 
production of ozone', ii) their atmospheric processing and subsequent deposition which leads to acidification and 
eutrophication in ecosystems through the formation of nitric acid (HNOs)'''^, iii) their role as precursors of the 
nitrate radical (NO3) (which controls the oxidizing capacity of the nocturnal atmosphere'"''), and iv) their effect on 
the radiative forcing of the atmosphere, directly through the absorption of radiation by NO2, or indirectly through 
ozone formation and altering the lifetimes of several greenhouse gases*''. Although NOx ^re also naturally 
produced by soil emissions, biomass burning and lighting, their main source is anthropogenic through the 
combustion of fossil fuels such as coal, oil (mainly in urban traffic) and gas. Hence, NOx emissions in a country 
or region are directly related to its industrial and economic development, and are consequently affected by 
changes in the economic and industrial activity, in addition to the implementation of environmental policies'""'^. 

Most of the NOx emissions from combustion systems are in the form of NO, although in diesel vehicles NO2 
emissions can reach up to 70% of the emitted NOx'"*. From 1996 to 201 1, the gasoline/gas-oU percentage ratio in 
Spain has changed from 75/25 to 45/54(1% corresponds to electric vehicles)''', therefore a relatively higher 
amount of NOx emissions is currently produced in the form of NO2. As sunlight produces an equilibrium between 
NO and NO2 through a photochemical cycle involving ozone', space-based observations of NO2 can be used to 
locate sources of NOx emissions due to the short lifetime of N02'^. The first studies about the temporal evolution 
of tropospheric NO2 were reported by Richter et al. for the 1 996-2004 time period in the main industrialized and 
populated areas of the world, and Irie at al." for 1996-2002 focusing on East Asia. The first study reporting the 
effect of the economic downturn in the NO2 levels was published by Lin and McElroy'^, showing a clear decrease 
in China from late 2008 to late 2009. Since these studies there have been other reports about the evolution of 
tropospheric NO2 as seen from space'" More recently, the global study on tropospheric NO2 trends by Richter 
et al."" has been updated for the 1996-2011 period by HUboU et al.''\ Previous reports of NO2 trends have also 
highlighted a clear reduction in the NO2 levels over urban areas due to both the application of anthropogenic NOx 
emission control policies and the effects of the global economic crisis'"'"'^'. 
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In this work we report on the NO2 geographical distribution and 
concentration trends over Spain using combined multiple satellite 
datasets and ground-based observations from five metropolitan 
areas in the country: Madrid (3.233.527 inhabitants), Barcelona 
(1.620.943 inhabitants), Valencia (797.028 inhabitants), Sevilla 
(702.355 inhabitants) and Bilbao (351.629 inhabitants)". Ground- 
based observations in Leon and Ponferrada are also reported, as these 
relatively small cities are located in the region where satellite mea- 
surements detect the higher decrease. We assess the contribution of 
environmental policies and the recent severe economic recession to 
the observed reduction in NO2 levels over Spain in the last two 
decades. Finally, the variation of the main economic indicators of 
the country is compared with the measured NO2 levels, to infer the 
possible relationship between the NOx emissions and the economic 
evolution of Spain. 

Results 

Satellite-based NO2 observations and trends. Figure 1 shows averaged 
tropospheric NO2 Vertical Column Density (VCD) maps for 1996 
and 201 1 over the Iberian Peninsula, with data from the GOME and 
SCIAMACHY instruments, respectively. Due to the different ground 
pixel size between these different sensors, a resolution correction 
method is applied to the low spatial resolution data from GOME'"" 
(see Methods). This procedure allows building a time-consistent 
NO2 dataset from 1996 to 2012. The figure includes annual, winter 
(Dec-Jan-Feb), and summer (Jun-Jul-Aug) NO2 column averages for 
the first and last years of the time series. Data from the GOME 
resolution corrected data product have been used for 1996. The 
plot for the last year of the series is made from 2011 SCIAMACHY 
measurements, as the operative cycle of this instrument ended in 
April 2012. Figure 1 clearly illustrates the difference between 
concentrations and geographical distribution of tropospheric NO2 
in 1996 and 2011. The seasonality of NO2 values is also evident for 
both years, showing higher concentrations in winter that contrast 
with summer values. As expected in mid-latitude developed 
countries, the lower solar radiation levels and the higher energy 
consumption, due to building heating systems, produce higher 
values of tropospheric NO2 in the winter season"^'. 

The highest levels of NO2 are observed in the regions of Madrid 
and Cataluiia, the two largest urban areas of Spain. Madrid's eco- 
nomy is based on the service, construction and industry sectors. Its 
location in the centre of the country also makes Madrid the main 
transport knot within the Iberian Peninsula. Cataluiia is located in 
the north-east coast of Spain, and its economy is mainly based on 
service, tourism and industry. Most of the economic activity in these 
two regions is located close to the metropolitan areas of the cities of 
Barcelona and Madrid. 

The north and east coastal regions of the peninsula are also areas 
with considerable industrial activity, and is therefore reflected in high 
NO2 levels (Fig. 1). The north-west part of the peninsula (Asturias, 
Leon and Palencia provinces, see Figure SI for their location) is an 
important region of mining and power generation industries. Also, 
considerable maritime transport, mining and steel industries activ- 
ities take place in the north coast area of the Gulf of Vizcaya, close to 
the city of Bilbao. At the east coast, most of the economic activity is 
located in the Valencia region and the Ebro basin, whereas the 
Guadalquivir basin concentrates most of the economic and industrial 
activity in the south-west part of Spain, where cities such as Cordoba, 
Sevilla and Huelva are located. 

The reduction in the NO2 levels over recent times is evident when 
comparing the first and last years of the satellite dataset (Fig. 1). The 
decrease is clear in the annual and summer averages, although the 
change in winter is less pronounced. Both annual and summer 
averages show a marked reduction in Madrid, Cataluiia, and the 
north-west part of the peninsula. In winter the drop is less apparent, 
but still appreciable in Madrid, east coast and north-west areas. Note 



that the 2011 winter was particularly cold'^* and thus could be 
responsible for a higher use of heating systems, incrementing NOx 
emissions during that winter. The more widespread NO2 distribu- 
tions over the regions of Madrid and Cataluiia in 1996 are reduced 
and limited by 20 1 1 to the metropolitan and localised industrial areas 
around the cities of Madrid and Barcelona. The decrease in NO2 
levels and its geographical spread is also clear along the north and 
east coasts, and Ebro and Guadalquivir basins (Fig. 1). 

To study the time evolution of NO2 levels over the Iberian pen- 
insula, the trend model described in HilboU et al.'^ has been applied 
here to the NO2 measurements from the GOME and SCIAMACHY 
instruments. Figure 2 shows the integrated evolution of tropospheric 
NO2 levels during the 1996-2011 time period, including absolute and 
relative changes. Both plots illustrate the decrease in NO2 concentra- 
tions over a large part of the peninsula, and particularly in the main 
economic and industrial areas. NO2 levels in the area of Madrid have 
decreased on an annual average by up to 1.9 X 10" molecules cm""^, 
which adds up to 53% of the levels measured in 1996. This reduction 
is much higher than the 20% previously reported for a different 
period (2004-2010) using OMI measurements by CasteUanos and 
Boersma, 2012'°. This difference is due to the different time period 
considered and also note the difference in local time at which differ- 
ent instruments overpass Spain. The measurement times for GOME 
and SCIAMACHY are -10:30 and -10:00 local time, while OMI is 
-13:45. GOME and SCIAMACHY measurements are closer to the 
diurnal NO2 peak in the early morning, when NOx emissions are 
more affected by transport sources at the rush hour. Emissions from 
industry and power plants have a less pronounced diurnal variation; 
consequently lower NO2 VCDs are expected from the OMI measure- 
ments resulting in lower relative decreases. These differences in the 
SCIAMACHY product compared to OMI have been previously 
studied^'"^', reporting higher NO2 levels from SCIAMACHY in the 
northern midlatitudes and regions where the NOx emissions come 
mainly from fossil fuel sources. Other reasons are the different pixel 
size of OMI and the differences in the instruments and retrieval 
algorithm. Considerable reductions of —40, —36, and —51% are 
measured in the areas of Barcelona, Bilbao and Valencia, respect- 
ively. In Sevilla, however, a smaller decrease of 27% is observed. 

The most widespread and remarkable reduction in NO2 levels is 
detected in the northwest industrial area and the north coast, with the 
highest decreases in the Asturias region and the provinces of Leon 
and Palencia (Fig. 2). These regions have experienced a large reduc- 
tion in the industrial, mining and coal power plants activities, which 
together with the application of new emission control technologies 
have resulted in a drastic reduction in NO2 concentrations. The 
decline in mining and industrial activities in these regions from 
2004 to 2009, have led to a fall in the NOx emissions of -52% from 
the main power plant stations of the region""*. A summary of the 
satellite NO2 trends by geographical areas is given in Table 1, includ- 
ing absolute and relative changes. 

Recently, CasteUanos and Boersma'" for Europe and Vrekoussis 
et al."^' for Greece, have also reported a general negative trend in NO2 
concentrations accelerated since the beginning of the economic 
recession (i.e., 2008). The trend model can also be applied to the 
pre-recession and recession time periods, although the nonlinearities 
in the general trend and the short length of this individual periods 
yield few significant data. Nevertheless the purpose of this work is to 
show the evolution of the NO2 levels rather than specifically analyse 
the short time period trends. Overall, from 1996 to 2011, the inte- 
grated NO2 VCD as observed from space over the entire Spanish 
territory has been reduced by —22%. 

NO2 trends as measured by in-situ air quality networks. To quan- 
tify the variation of surface NO2 levels in the cities of Barcelona, 
Bilbao, Madrid, Sevilla and Valencia, data from their ground air 
quality monitoring networks (AQMN) are used. All the in-situ 
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Figure 1 | Averaged tropospheric NO2 VCD for different seasons during 1996, using the GOME resolution-corrected dataset (a, c and e), and 2011/ 
2012 using SCIAMACHY data (b, d, f). Top: annual averages. Middle: summer (Jun-Jul-Aug) averages. Bottom: winter (Dec-Jan-Feb) averages. 
Figure created by the authors using Matlab R213a. 



monitoring stations employ the chemiluminiscence method 
described in the Methods section. Figure 3 shows the annually- 
averaged NO2 levels from 1996 to 2012 measured by the available 
AQMN. All cities show a general decrease in the NO2 levels during 
the last two decades. Regression analysis yields considerable absolute 



reductions of 67% for Valencia and 52% in SeviUa, while Barcelona 
exhibits a shallower 18% decrease. Madrid and Bilbao show 
intermediate decreases of 37% and 35%, respectively. Overall, the 
average of all cities considered here shows a decrease of 44% in the 
NO2 concentrations. 
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Figure 2 | Change in the NO2 tropospheric VCD over the Spanish territory derived from the shift trend model for the years 1996-2011 (data with 
significant trends at the 95% level). Figure created by the authors using Matlab R213a. 



The results of the statistical regression analysis performed for 
different cities and time periods are shown in Table 2, including 
absolute and relative changes. This summary also includes the com- 
bined average of all studied cities. Barcelona, Madrid, and Valencia 
show two well differentiated slopes for the pre-recession and reces- 
sion periods ( — 0.5, —0.7, and — 3.1/|ig m""* yr"' in Barcelona, 
Madrid and Valencia, respectively, during the pre-recession period, 
and —0.9, —4.3, and — 3.9/|.ig m"' yr"', respectively, during the 
recession period). In these three cases the slopes are higher during 
the recession period. These slopes correspond to annual decreases of 
0.9%, 1.1% and 4.3 for Barcelona, Madrid and Valencia, respectively, 
during the pre-recession period, and 1.9%, 7.8% and 9.1%, respect- 
ively, during the economic recession. SevUla shows similar slopes for 
both periods, being slightly higher during the pre-recession period 
(— 1.7/|.ig m"^ yr"') than during the economic recession ( — 1.5/|ig 
m"'' yr"'). By contrast with the other cities, Bilbao shows a higher 
slope during the pre-recession period than during the economic 
recession ( — 1.4 and — 0.8/|j,gm"'yr"' respectively). Overall, all cities 
except BUbao, show higher reductions during the recession period 
than during the pre-recession. This result suggests an initial reduc- 
tion driven by the progressive implementation of environmental 
policies and technology, whereas the second slope can be attributed 
to the rapid decrease of the economic and industrial activity during 
the economic downfall between 2008 and 2012. 

In order to confirm the pronounced decline in NO2 levels reported 
by the satellite measurements in the North-West region, AQMN data 
in the cities of Ponferrada and Leon (province of Leon) were also 



investigated. These data show a decrease of 80% in Ponferrada and 
66% in Leon since 1997. Both cities also show a higher decrease 
during the economic recession period, 72% and 75%, compared to 
63% and 41% during the pre-recession period for Ponferrada and 
Leon, respectively. 

The general annually- averaged NO2 trend from both satellite and 
ground stations is shown in Figure 4. The satellite averages have been 
made integrating over the entire Spanish territory in the Iberian 
Peninsula, while the ground stations data correspond to yearly aver- 
age in the studied cities. Hence the smoother shape in the AQMN 
dataset, which is more affected by local effects than the satellite 
measurements. Note that both datasets are not directly comparable, 
although their combination provides a clear picture about the evolu- 
tion of NO2 levels in Spain during the last two decades. 

Both datasets show two distinct slopes corresponding to the pre- 
recession and recession periods, the slope for the last period being 
higher. The purpose of this work is to show how the NO2 levels have 
changed in absolute terms during the last two decades, rather than to 
perform a statistical trend study. Therefore we refer here to 'general 
trends', despite the evident variable nonlinear trends observed in the 
time series. The nonlinearity in the NOx emission trends have 
already been reported by Konovalov et al.^' for the case of Madrid. 

Discussion 

In order to gain insights about how the economic activity affects NOx 
emissions in Spain, the levels of NO2 are compared to the main 
economic indicators, including Gross Domestic Product (GDP), 



Table 1 | Absolute and relative changes in the NO2 tropospheric VCD from 1 996-201 1 , derived from the trend model for different 
geographic regions. (See Table SI for latitude-longitude coordinates of the studied regions) 



Geographic region 


Absolute significant 
changes/molecule cm"^ 


Error/ ±mol 


ecule cm ^ 


Relative significant 
changes/% 


Error/±% 


East coast of Valencia 


-1,27E + 15 


3,60E 


+ 14 


-33 


9 


Ebro Basin and Cataluiia 


-1,34E + 15 


4,97E 


+ 14 


-23 


9 


Area of Madrid 


-3,10E + 15 


9,48E 


+ 14 


-53 


16 


Area of Sevilla 


-1,16E + 15 


5,33E 


+ 14 


-27 


12 


Area of Barcelona 


-2,98E + 15 


9,89E 


+ 14 


-40 


13 


Area of Bilbao 


-2,25E + 15 


7,29E 


+ 14 


-36 


12 


North-West 


-2,87E + 15 


9,22E 


+ 14 


-46 


15 


Gulf of Vizcaya 


-1,34E + 15 


6,24E 


+ 14 


-24 


11 


Guadalquivir basin 


-6,36E + 14 


3,86E 


+ 14 


-17 


10 


Area of Valencia 


-2,21E -H 15 


5,89E 


+ 14 


-51 


13 


Integrated over Spain 


-9,47E + 14 


3,55E 


+ 14 


-22 


8 
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Figure 3 | NO2 annual averages from air quality monitoring stations. 



Industrial Production Index (IPl, base 2005) and Oil Consumption. 
The GDP is the market value of all officially recognized final goods 
and services produced within a country in a given period of time. The 
IPI measures the real production output of manufacturing, mining 
and utilities of a country. The Oil Consumption refers to the number 
of barrels consumed in a country per day. NOx natural emissions 
account for 2.8-3.8% of total emissions in Spain (National Emission 
Inventory™) during the studied period. Therefore the trends in eco- 
nomy and industrial activities can be considered the main source of 
variation for the NOx emissions in Spain. Bottom-up emissions data 
from the National Emission Inventory"*" and the Centre on Emissions 
Inventories and Projections" are shown in Figures 82 and S3. Note 
however that the trends in these emission inventories, which foUow 
the trends in the economic parameters of the country, show a dis- 
crepancy with both on-ground and satellite data (see Figs. 5, S2 and 
S3). This highlights potential uncertainties in the trends of emission 
inventories and suggests that further research is needed to reconcile 
differences in the trends of observed NO2 and the emission invent- 
ories during the pre-recession period in Spain. 

Although Spain has failed to achieve the Gothenburg protocol 
targets, it has reduced the NOx emissions by about 20% from 1990 
to 2011, as reported by the European Environment Agency'"'. This 
fact is reflected in the NO2 measurements through a reduction by 
22% in the NO2 levels integrated over Spain as observed from space, 
and by an average of 45% in the AQMN observations in the five cities 
studied here. This reduction is considerably higher in the main 
industrial and economic cores of the country, such as Madrid and 
the north-west part of the Iberian Peninsula, with reductions of up to 
50% in the areas of Madrid and Valencia. In general, there is a 
correspondence between the tendencies of on-ground and satellite 



data (Figure 4), having two well differentiated slopes in the trends 
(i.e. pre-recession and economic recession periods). 

On the other hand, the main economic indicators were gradually 
increasing since 1996 until the beginning of the economic crisis in 
2008, when all indicators suddenly fell. These economic metrics refer 
to the whole country, so satellite NO2 VCDs integrated over the 
whole Spanish territory in the Iberian Peninsula are used to assess 
the correlation between economic fall and NO2 concentrations. In 
general we find a good correspondence between the economic indi- 
cators and the evolution of NO2, as shown in Figure 5, despite the 
2001-2007 period in which the trends are different. The NO2 
increase from 1996 to 2000 coincides with a period of economic 
growth marked by a huge increase of the IPI and a small, but appre- 
ciable, increase in the GDP. During this first period it can be con- 
sidered that the economic rise is directly reflected in the NO2 levels. 
From 2000 onwards the IPI slows down its rising while the NO2 
VCDs faU from 4.3 X 10'= to 3.6 X 10'= molecules cm~\ and the 
GDP undergoes a huge and constant growth until the start of 
the recession in 2008. Consequently, during the period in which 
the economic activity of the country is rising, the application of 
environmental policies and the optimization in the combustion pro- 
cesses are likely responsible for the reduction in the emissions of NOx 
measured from 2000 onwards. FinaUy during the recession (2008- 
2012) both indicators and the NO2 levels decrease considerably, 
showing a clear effect of the economy on the levels of anthropogenic 
pollutants. 

The same information can be extracted from the oil consumption 
evolution in Spain (Figure 5). From 1996 to 2008 it has been growing 
at constant rate, starting a deep decrease in 2008 with the recession. 
Note that the 33% increase in the oil consumption from 1996 to 2008 
should have led to an increase in the NO2 levels during this time 
period. This indicates the positive effect of emissions legislation upon 
combustion technologies. 

Our results show a lack of correlation between NO2 levels and 
economic indicators from 2001 to 2007, however Vrekousis et al.^' 
reported good agreement between NO2 and economic markers in 
Greece from 2000 to 2011. This difference may be due to the fact that 
Greece has not reduced its NOx emissions, as a result of envir- 
onmental legislation, at the same rate as Spain '\ and hence the better 
correspondence between economic evolution and NOx emissions in 
Greece. This reduction in the NO2 levels in Spain during the 2001- 
2007 period, when all the economic factors increase, is a clear signal 
that the application of environmental policies are efficiently reducing 
the NOx emissions. This has also been reported for other cities of the 
European Union as shown by Castellanos and Boersma'". 

To summarize, the datasets presented here can be divided into two 
differentiated parts: i) from 1996 to 2008 the application of envir- 
onmental policies reduced the NOx emissions, with a slight increase 
from 1996 until 2000 before the effects of environmental policies 
start to pay off; ii) during the 2008-2012 recession period the eco- 



Table 2 | Regression analysis of NO2 trends as measured by air quality monitoring stations. Slopes 


in Table S2 








1996-2012 


1 


996-2008 


2008-2012 


City 


Decrease/% 


yr"' Absolute decrease/% 


Decrease/% yr 


' Absolute decrease/% Decrease %yr ' 


Absolute decrease/% 


Barcelona 


-1.16* 


-18.60 


-0.92 


-1 1.99 


-1.89* 


-7.56 


Bilbao 


-2.19* 


-35.02 


-2.89 


-37.59 


-2.27* 


-9.08 


Madrid 


-2.20 


-37.45 


-1.1 1 


-14.42 


-7.84 


-39.18 


Sevilla 


-3.28** -52.43 


-3.64** 


-43.66 


-4.77 


-23.87 


Valencia 


-4.47 


-67.09 


-4.32 


-47.50 


-9.13 


-45.66 


Average*** 


-2.64 


-44.93 


-2.23 


-29.02 


-6.86 


-34.29 


*No data available 


or 201 2. 












**No data available 


for 1996. 












***Average of the five cities 
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nomic crisis remarkably reduced the emission even more than the 
environmental policies. These trends are observable through both 
datasets, although the absolute changes are higher in the ground- 
based measurements, showing decreases up to 67% in the city of 
Valencia. These differences between the two datasets are, as 
explained in the Methods section, due to the different type of mea- 
surements considered. The ground stations represent surface obser- 
vations within the cities, while the satellite measurements result from 
integrating the tropospheric column over the whole ground pixel of 
the instrument, which is larger than the largest of the studied cities. 
Thus we regard the ground measurements as representative of the 
local NO2 levels and trends in the cities, while the satellite-based 
observations represent column-integrated NO2 levels over larger 
areas around the cities. The 20% emissions reduction reported by 
the European Union agree better with the satellite dataset, as these 
emissions reductions refers to the whole country rather than to spe- 
cific cities or regions. 



The lack of correlation between the economic parameters of Spain 
and the levels of NO2 during the 2001-2007 period, with NO2 levels 
decreasing as the GDP, IP and Oil Consumption grow until 2008, 
shows that the improvement in the combustion technologies and 
environmental policies implemented to reduce the nitrogen oxides 
emissions are paying off Therefore the air quality is improving, 
although there is still margin for further improvement as demon- 
strated by the fact that the daily thresholds of NO2 and O3 levels are 
often exceeded in Spanish cities. The steeper slope in the NO2 reduc- 
tion during the recession, when the economy falls, is clear in both 
datasets. Hence for the case of Spain, still on recession during 2012, it 
can be expected that the slope of the decrease will recover to the pre- 
recession values once the economy starts growing again. Finally, this 
decline in the NO2 levels during the last two decades has led to an 
overall improvement of the air quality in the Iberian Peninsula, and is 
likely affecting the chemistry of the troposphere in three possible 
ways: 1) with the reduction in the NOx levels, a variation in the 
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Figure 5 | Satellite tropospheric NO2 VCD averaged over the Spanish territory, along with the main economic metrics of Spain, Gross Domestic 
Product (GDP), Industrial Production Index (IPI) and Oil Consumption. 
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tropospheric ozone levels is also expected, 2) this denoxification of 
the atmosphere has probably led to lower nitrogen inputs to ecosys- 
tems, and 3) the levels and seasonality of HNO3 formation as influ- 
enced by traffic emissions^^^ have likely been affected. 

Finally, it is worth mentioning that future monitoring of air pol- 
lution in megacities and industrial regions will greatly benefit from 
geostationary satellite missions, which will considerably enhance the 
temporal and spatial resolution of the observations. Therefore it will 
provide a more accurate tool to verify the effect of implementing 
environmental policies and to follow the atmospheric impacts of 
changing emissions. 

Methods 

Description of datasets. Two main datasets including space-based observations and 
air quality monitoring networks (AQMN) measurements are used to study the 
evolution of tropospheric NO2 over Spain. The satellite observations provide an 
overview of the NO2 geographical distribution and variation not only over major 
urban areas, but also over the entire Iberian Peninsula, while in-situ ground 
monitoring stations provide insightful information of the pollutant levels evolution 
within cities. Note that the spatial resolution of current satellite -based spectrometers 
is still too low to properly assess detailed NOx emissions at a sub-urban scale or from 
medium-size cities. On the other hand AQMN provide information of pollutants at 
surface level, while satellite based measurements integrate the NO2 vertical column 
throughout the troposphere (once the stratospheric contribution has been removed). 
Therefore, both types of measurements are not directly comparable but 
complementary. 

The satellite dataset is based on measurements from the GOME and SCIAMACHY 
instruments, using the trend model developed by HilboU et al.' \ This model accounts 
for the ground pixel size and instrument differences between GOME and 
SCIAMACHY, as described below in this section. The in-situ dataset has been built 
using data from the AQMN in each of the cities included in this exercise. Data from a 
total number of --55 monitoring stations of urban and sub-urban type have been used 
during the period of study, with different typology and locations representative of the 
major emission points within each individual city. These AQMN are operated and 
maintained by the regional governments of the country, according to current 
European legislation. 

Surface measurements. The method used to detect NO2 in the monitoring station is 
described in EN 14211:2005 'Ambient air quality - Standard method for the 
measurement of the concentration of nitrogen dioxide and nitrogen monoxide by 
chemiluminescence, defined by the European Directive (2008/50/CE) as reference 
method. In this method a pre-reactor separates the NO-NO2 reaction from the 
background chemiluminescence, allowing accurate auto-zeroing of the analyzer. 
The equipment works under a rigid maintenance program being periodically tested 
and calibrated. 

Satellite measurements. In this work measurements performed by the satellite based 
instruments GOME^' on board ERS-2, and SCIAMACHY-''^ " on board ENVISAT are 
employed. The ground pixel of these instruments is 40 X 320 km^ and 30 X 60 km^ 
for GOME and SCIAMACHY, respectively. The raw data have been analysed as 
described in Hilboll et al., 2013a''^. In short, differential optical absorption 
spectroscopy'** in the 425-450 nm wavelength window'^ is applied to the measured 
radiances. The influence of stratospheric NO2 has been corrected for using the 
algorithm detailed in Hilboll et al. 2013"^°, using stratospheric NO2 fields from the 
BBdCTM model^'"^\ Tropospheric air mass factors have been calculated with 
the radiative transfer model SCIATRAN'*^. 

Correction for the difference in ground pixel size between GOME and 
SCIAMACHY. To show the spatial distribution of tropospheric NO2 columns (see 
Fig. 1) from GOME measurements, the retrieved data have been post-processed to 
account for the very low spatial resolution of the GOME instrument. The method is 
based on spatial averaging of earthshine spectra measured by SCIAMACHY and the 
extraction of a spatial pattern of the effect of ground pixel size. This effect of ground 
pixel size is subsequently super-imposed on the original GOME measurements, 
yielding a data set "GOME res.corr". In this data set, the magnitude of the 
tropospheric NO2 columns is retrieved from GOME measurements, while the spatial 
structure is derived from the full SCIAMACHY data set 2003-2011. The method is 
described in Section 4 of Hilboll et al., 2013^^ 

Trend model accounting for the difference between GOME and SCIAMACHY 
measurements. The rate of linear change of tropospheric NO2 columns has been 
calculated using the algorithm described in Section 5 of Hilboll et al., 2013^'\ In short, 
the trend model explicitly accounts for (a) an offset between the two instruments and 
(b) for a change in the amplitude of the seasonal variation. The time series of monthly 
averages of NO2 measurements Y (t) is described by the equation 

Y{t)^^i-\-cot^cWit)-\-Sit)-^N{t) (1) 
where p is the VCDti-op measurement at time t — 0, w is the monthly trend 



component, and t is the time in months since January 1996. S is the level shift between 
GOME and SCIAMACHY measurements occurring at time t — Tq (which we set to 
January 2003), and U(t) is the step function 

UitU'- (2) 
I 1, t>To ^ ^ 

The seasonal component S(t) has the same shape (represented by harmonic 
functions) and varying amplitude for the two instruments, and is modelled by four 
sine and cosine components, each. The noise process N(t) is assumed to be 
autoregressive with lag 1, and only trends which are significant at a 95% level are 
considered. 
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